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Abstract: In forest ecosystems, interactions between overstory trees and understory herbs play an 
important role in driving plant species diversity. However, reported links between overstory tree and 
understory herb species diversity have been inconsistent, due to variations in forest types and 
environmental conditions. Here, we measured species richness (SR) and diversity (Shannon- Wiener (H') 
and Simpson's (D) indices) of overstory trees and understory herbs in the protected Tianshan Wild Fruit 
Forest (TWFF), Northwest China, to explore their relationships along the latitudinal, longitudinal, 
elevational, and climatic (current climate and paleoclimate) gradients in 2018. We found that SR, and H' 
and D diversity indices of overstory trees and understory herbs exhibited a unimodal pattern with 
increasing latitude and elevation (P<0.05) and negative associations with longitude (P<0.01). Along the 
climatic gradients, there were U-shaped patterns in SR, and H' and D diversity indices between trees and 
herbs (P<0.05). SR, and H' and D diversity indices for overstory tree species were positively associated 
with those for understory herbs (P<0.01). These findings indicate that overstory trees and understory 
herbs should be protected concurrently in the TWFF to increase effectiveness of species diversity 
conservation programs. 
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1 Introduction 


Diversity in forest ecosystems reflects interactions between overstory trees 
and understory herbs (hereafter referred to as trees and herbs, respectively) 
(Gilliam, 2007; Pan et al., 2013) and is regulated by a series of feedbacks 
between the two strata. For example, trees determine light availability and 
heterogeneity of soil fertility available to herbs (Muller, 2003; Neufeld and 
Young, 2003), while herbs drive positive and negative feedback processes with 
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trees, such as competition for resources among herb and tree seedlings and 
the return of nutrients to the soil from herb litter (Yarie, 1980; Gilliam, 2007). 
Thus, it is important to understand the relationship between tree and herb 
species diversity in forests to improve practical management and conservation 
of forest biodiversity. 

The relationship between tree and herb species diversity has been studied in 
a range of forest types, such as temperate deciduous, subalpine forest, and 
mixed conifer forests (Ewald and Freising-Weihenstephan, 2002; Ingerpuu et 
al., 2003; Gilliam, 2007; Vockenhuber et al., 2011), and was found to vary 
among the types. In global temperate and German deciduous forests, tree 
species diversity was positively associated with that of herbs (Gilliam, 2007; 
Molder et al., 2008; Vockenhuber et al., 2011). While in mountain forests of 
Bavarian Alps and deciduous forests of North America, there was no 
relationship between tree and herb species diversity (Ewald and Freising- 
Weihenstephan, 2002; Houle, 2007). Besides, high levels of tree species 
diversity in conifer forests were negatively associated with diversity of herbs, 
due to low levels of soil pH caused by conifer species (Berger and Puettmann, 
2000). 

This inconsistency in the relationship between tree and herb species 
diversity among forest types may be result of a number of factors, including 
stand age and diversity condition, and environmental gradients. For example, 
the relationship between tree and herb species diversity is known to vary with 
stand age in hardwood forests, where there is an association in mature stands, 
but not in young stands (Gilliam et al., 1995); herb diversity tends to be greater 
in stands of mixed rather than single species (Barbier et al., 2008); and the 
relationship between tree and herb species diversity varies along the 
environmental gradients (Gilliam and Robert, 2003; Vockenhuber et al., 2011). 
The effect of environmental gradients on forest plant diversity is complex. For 
example, at the global scale, tree and herb species diversity is unimodal with 
increasing elevation, but their diversity peaks appeared at different elevations 
due to differences in physiological tolerance and niche partitioning (Guo et al., 
2013). Also, it is reported that in the eastern Himalayas, herb species diversty 
is positively related to annual levels of precipitation, and there is a unimodal 
relationship with increasing temperature, whereas the relationships for tree 
diversity are the reverse (Kluge et al., 2017). Differences in evolutionary history 
and sensitivity to local climate have led to contrasting tree and herb species 
diversity responses to latitudinal gradient in forests in China (Wang et al., 
2012), indicating that the relationship between tree and herb species diversity 
is environment context-dependent. Thus, further research is needed to clarify 
environmental gradient effects on tree and herb species diversity relationships. 

The Tianshan Wild Fruit Forest (TWFF), which originated from the late 
Tertiary and early Quaternary, is a special broad-leaved forest and a priority 
conservation ecosystem in China, comprising a mixture of paleotemperate 
broad-leaved forests and north forest meadows (Chang, 1973; Xu et al., 2006). 
A total of 441 vascular plant species have been identified in the TWFF, 
including many that are endangered and relic plants (Yang et al., 2003), such 
as Malus sieversii and Prunus sogdiana that are wild ancestors of many 
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cultivated tree fruit species (Yang et al., 2003; Duan et al., 2017). Although 
plant diversity in the TWFF has been widely reported (Yang et al., 2003; Li et 
al., 2011; Fang et al., 2019), however, our understanding about variation of 
diversity along the environmental gradients and links between tree and herb 
species diversity remains limited, as these studies neglecting effects of 
paleoclimate changes on diversity in the TWFF and pooling tree and herb 
species together as total forest diversity (Yang et al., 2003; Fang et al., 2019). 
Therefore, in this study, we surveyed tree and herb species in the TWFF to 
quantify changes in diversity with latitude, longitude, elevation, and climate, 
and to clarify the relationship between tree and herb diversity in this forest 


type. 
2 Materials and methods 


2.1 Study area 


The study was conducted at five sites with contrasting tree and herb species 
composition and environmental factors in the lli Kazak Autonomous Prefecture, 
Xinjiang Uygur Autonomous Region in Northwest China. The study area is 
located between Tianshan Mountains and Altay Mountains (Fig. 1), which is 
the center of the TWFF in China (Chang, 1973; Yang et al., 2003). The 
Xinyuan and Gongliu study sites are located at the southern edge of Tianshan 
Mountains and are similar in tree species composition (Fig. 1), with pure 
stands of Malus sieversii, but the two study sites differed in herb composition, 
where Urtica dioica, Festuca gigantean, and Dactylis glomerata dominated at 
the Xinyuan study site and Cannabis sativa, Poa annua, and Setaria viridis 
dominated at the Gongliu study site. Study sites of Huocheng, Tuoli, and Emin, 
which are located at the northern edge of Tianshan Mountains and southern 
edge of Altay Mountains (Fig. 1), are comprised by mixed stands, dominated 
by Malus sieversii, with abundant Juglans regia, Armeniaca vulgaris, Prunus 
sogdiana, Crataegus chlorocarpa, and Sorbus tianshanica, while dominant 
herbs comprised Impatiens brachycentra, Achillea millefolium, Carex 
pediformis, Potentilla bifurca, Bromus japonicas, and Iris tenuifolia, depending 
on site condition. The study area ranges from 843.8 to 1406.0 m a.s.l. in 
elevation and is characterized by a temperate continental climate, with mean 
annual precipitation (MAP) ranging from 196 to 266 mm and mean annual 
temperature (MAT) ranging from 3.2°C to 7.0°C (Table 1); mountain black 
brown soils dominate the area (Lin and Cui, 2000). 
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Fig. 1 Distribution of the study sites 


Table 1 Geographic location, elevation, mean annual temperature (MAT) and mean annual precipitation 
(MAP) of the study sites 


Sites Latitude (°N) Longitude (°E) Elevation (m) MAT (°C) MAP (mm) 
Emin 46.37 83.99 1255.0 3.2 266 
Tuoli 46.15 83.55 843.8 5.7 235 
Huocheng 44.43 80.79 1173.2 7.0 196 
Gongliu 43.24 82.72 1261.0 6.4 196 
Xinyuan 43.38 83.60 1406.0 5.6 225 


2.2 Plant diversity 


Tree and herb species diversity was assessed in August 2018, when the herb 
diversity reached the maximum level. To measure diversity of tree species, we 
randomly placed five 20 mx20 m plots at each study site, with the exception of 
Tuoli where four plots were established. The 20 mx20 m sample size was 
recommended and used for monitoring networks of forest biodiversity in China 
(Mi et al., 2016; Wang et al., 2018). Tree and shrub species abundance was 
recorded in the plots; shrubs were classified as trees because shrubs in this 
area are rare and have a high height. We recorded herb species abundance 
within each plot in three (at the study sites with low diversity) or five (at the 
study sites with high diversity) randomly placed 1 mx1 m quadrats. We 
calculated species richness (SR), and Shannon-Wiener (H') and Simpson's (D) 
diversity indices for the recorded tree and herb species (Mölder et al., 2008; 
Vockenhuber et al., 2011), and classified herbs to functional group (forbs or 
graminoids) to analyze their relationships with tree species diversity 
(Vockenhuber et al., 2011). 


2.3 Environmental variables 


In each study plot, geographic information (longitude and latitude) and 
elevation were recorded using a portable GPS (eTrex H, Taiwan, China). We 
also extracted current climate variables (MAP, precipitation of coldest quarter 
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(PCQ), MAT, and mean temperature of coldest quarter (MTCQ)) and 
equivalent paleoclimate change anomalies since the Last Glacial Maximum 
(LGM) (MAPano, MATano, MTCQano, and PCQano, which are anomaly 
values of MAP, MAT, MTCQ, and PCQ since the LGM, respectively), as they 
affect trees and herb species diversity (Currie and Paquin, 1987; Feng et al., 
2016) and current diversity patterns, respectively (Jansson, 2003; LU et al., 
2018). Current climate variables were extracted from WordClim (1970—2000; 
http://www.worldclim.org/) at a resolution of 1 kmx1 km (Hijmans et al., 2005), 
and MAPano, PCQano, MATano, and MTCQano were calculated by 
subtracting their values at the LGM, which were extracted from the Community 
Climate System Models v3 (http://www. worldclim.org/paleo-climate1) at a 
resolution of 1 kmx1 km (Jansson, 2003; Otto-Bliesner et al., 2006), from 
corresponding current climate values (Feng et al., 2016). 


2.4 Data analyses 


A general linear model (GLM) was adapted to explore trends in tree and herb 
species diversity along the environmental gradients (Ma et al., 2010). 
Environmental variables in this study were highly correlated with each other, 
that could lead to a miscalculation in the GLM due to multicollinearity; 
therefore, we used a principal component analysis (PCA) to produce a set of 
independent principal components (PCs) to represent the environmental 
gradients (Wang et al., 2012). The first principal component (PC1) reflected a 
tradeoff between temperature (MAT and MTCQ) and precipitation (MAP and 
MAPano) (Fig. 2). For example, there were low temperatures and high levels 
of precipitation at the Emin study site, whereas temperatures were high and 
levels of precipitation were low at the Huocheng study site (Table 1; Fig. 2). 
The second principal component (PC2) represented elevation-induced 
variations in paleoclimate and current climate (MTCQano, MATano, and PCQ; 
Fig. 2). For example, the Xinyuan study site was located at high elevation with 
high levels of MTCQano and low levels of MATano and PCQ, whereas the 
Tuoli 
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Fig. 2 Principal component analysis (PCA) of variation in environmental variables along the latitudinal 
gradient. MAP, mean annual precipitation; MAT, mean annual temperature; MTCQ, mean temperature of 
coldest quarter; PCQ, precipitation of coldest quarter. MAPano, MATano, MTCQano, and PCQano are 
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anomaly values of MAP, MAT, MTCQ, and PCQ since the Last Glacial Maximum, respectively. PC1 and 
PC2 represent the first and second principal components, respectively. 


study site was characterized by low elevation, with low levels of MTCQano and 
high levels of MATano and PCQ (Fig. 2). Together, PC1 and PC2 explained 
87.3% of total variation in environmental variables, so they were used in the 
GLM analysis of latitudinal and environmental trends in diversity (Wang et al., 
2012). We also used the GLM to explore the relationship between tree and 
herb species diversity, where both linear and quadratic relationships were fitted 
separately and the best model was selected using the Akaike Information 
Criterion (AIC) (Brunham and Anderson, 2002). All statistical analyses were 
performed in R software v2.15.1 (R Development Core Team, 2012). 


3 Results 


3.1 Trends in diversity along the latitudinal and longitudinal gradients 


We found that SR, and H' and D diversity indices of tree and herb species were 
unimodal along the latitudinal gradient (P<0.01; Figs. 3a—f), indicating that in 
the TWFF, the peak of diversity appeared at the middle level of latitude. 
However, SR, and H' and D diversity indices of tree and herb species were 
negatively associated with longitude (P<0.01; Figs. 3g—l), demonstrating that 
western sites had a high level of diversity while eastern sites had a low level of 
diversity. 

3.2 Trends in diversity along the elevational and climatic gradients 


SR, and H' and D diversity indices for tree and herb species showed U-shaped 
patterns with increasing temperature (MAT and MTCQ) and decreasing levels 
of precipitation (MAPano and MAP) (PC1 environmental gradient; P<0.001; 
Figs. 4a—f). In contrast, tree specie's SR, and H' and D diversity indices, and 
herb specie's SR and H' diversity index exhibited a unimodal pattern with 
increasing elevation and MTCQano, and decreasing MATano and PCQ (PC2 
environmental gradient; P<0.05; Figs. 4g—l). 


3.3 Relationship between tree and herb species diversity 


There was a consistent positive, linear relationship of richness (SR) and 
diversity (H' and D) between tree and herb species (P<0.01; Figs. 5a—c). When 
herbs were classified as forbs or graminoids, there was a positive relationship 
between tree and forb species diversity (P<0.01; Figs. 5d-f). However, there 
was no relationship between tree and graminoid diversity (P>0.05; Figs. 5g-i). 


4 Discussion 


4.1 Trends in tree and herb species diversity along the latitudinal and 
longitudinal gradients 


Diversity patterns along the latitudinal gradient is a key issue in ecology and 
biogeography (Hillebrand, 2004; Mannion et al., 2013), and studies have 
indicated they may be random or vary from unimodal, to positive or negative 
linear, depending on forest types and spatial scale (Gaston, 2000; Willig et al., 
2003). In the TWFF, richness (SR) and diversity (H' and D) of tree and herb 
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species with increasing latitude was unimodal (Figs. 3a—f), supporting previous 
studies of forests in China (Wang et al., 2011, 2012). This unimodal pattern of 
diversity with latitude was driven by temperature, particularly by the MTCQ and 
temperature seasonality that is the standard deviation of mean monthly 
temperature (Wang et al., 2011; Qian, 2013). Levels of MTCQ gradually 
decreased with increasing latitude, and it has been shown that high levels of 
MTCQ at lower latitudes reduce the proportion of temperate species, while low 
levels of MTCQ at higher latitudes decrease the proportion of tropical species 
that lead to low levels of species diversity at both lower and higher latitudes 
(Wang et al., 2011). The peak of species diversity at the mid-latitude was 
attributed to the coexistence of a large numbers of tropical and temperate 
species (Wang et al., 2011). This mechanism may also exist in the TWFF 
because MTCQ was significantly negatively correlated with latitude (P<0.001) 
and the unimodal pattern of tree and herb species diversity was also found 
with increasing MTCQano (Figs. 4g—). 

In the TWFF, we found that diversity of tree and herb species was negatively 
associated with longitude (Figs. 3g—l), supporting a recent study of forests in 
China (Wu et al., 2018). The high 
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Fig. 3 Trends in species richness (SR), and Shannon-Wiener (H') and Simpson's (D) diversity indices for 
tree and herb species along the latitudinal (a-f) and longitudinal (g—l) gradients. Estimated relationship 
(solid line) was fitted with confidence intervals (dashed lines) at P=0.05 level. 


levels of diversity at the western study sites and low levels at the eastern study 
sites may be attributed to variation in precipitation. For example, precipitation 
in Xinjiang in spring and summer 
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mainly derives from water vapor from the Mediterranean and Caspian Seas, 
and the Atlantic that 
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Fig. 4 Trends in SR, and H' and D diversity indices for tree and herb species along the PC1 (a-f) and 
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PC2 (g-l) environmental gradients. Estimated relationship (solid lines) was fitted with confidence intervals 
(dashed lines) at P=0.05 level. NS, no significant relationship. Positive (+) and negative (—) signs 
indicated change in direction of environmental variables along the PC1 and PC2 environmental 
gradients. 
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Fig. 5 Relationships of SR, and H' and D diversity indices of trees with herbs (a—c), forbs (d—f), and 
graminoids (g—i). Estimated relationship (solid line) was fitted with confidence intervals (dashed lines) at 
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gradually decreases eastwards (Yang and Liu, 2018). The negative 
relationship between tree and herb species diversity with increasing longitude 
in the TWFF may be partly driven by this unique precipitation pattern, as global 
diversity was also positively related with levels of annual precipitation (Gaston, 
2000). 


4.2 Trends in species diversity along the elevational gradient 


In the TWFF, tree species diversity was found to be unimodal with increasing 
elevation (Figs. 4g—l), supporting previous studies in Himalaya Mountains and 
temperate forests in Taibai, Taishan, and Laoshan mountains in China (Acharya et 
al., 2011; Zhang et al., 2016; Xu et al., 2019). However, the unimodal pattern of herb 
species diversity was in contrast to the negative and U-shaped pattern found in 
Taishan and Laoshan mountains (Zhang et al., 2016). A global study demonstrated that 
types of species diversity pattern along the elevational gradient were closely related to 
latitude and altitude, because unimodal patterns of species diversity have been 
recorded at high latitudes (>40°N) in the Northern Hemisphere (Guo et al., 2013). 
The TWFF lies between 43.23°N and 46.36°N, and the unimodal pattern of tree and 
herb species diversity was consistent with these global predictions. A previous study in 
the study area found a negative relationship between species diversity and elevation 
within the range from 1100.0 to 1500.0 m (Li et al., 2011). The greater range of 
elevation in our study (from 843.8 to 1406.0 m) increased the size of the dataset and 
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prediction power of the statistical model and likely led to the confirmation of a 
unimodal species diversity pattern (Guo et al., 2013). 

The unimodal pattern of species diversity along the elevational gradient was 
generally explained by environmental filtering and species competition 
(Lomolino, 2001). For example, low levels of species diversity at high 
elevations were mainly caused by decreases in temperature and availability of 
soil nutrients and increased isolation, because only few species may survive in 
such harsh environmental conditions (Guo and Berry, 1998; Lomolino, 2001). 
In contrast, low levels of species diversity at lower elevations may have 
resulted from highly competitive ability of native species (Guo and Berry, 1998; 
Lomolino, 2001). The peak in level of species diversity at mid-elevations can 
be attributed to suitable environmental conditions that allow for the greater 
survival and growth of a wider range of species (Guo and Berry, 1998; 
Lomolino, 2001). 


4.3 Trends in species diversity along the climatic gradients 


Paleoclimate, particularly paleoclimate change since the LGM, was found to 
elicit significant effects on contemporary patterns of tree and herb species 
diversity (Montoya et al., 2007; Feng et al., 2016). We found that diversity of 
tree and herb species in the TWFF initially decreased and then increased, with 
decreasing levels of MAPano and MAP (Figs. 4a-f), indicating that the climate 
has become more wetter and drier with greater levels of species diversity since 
the LGM. At the time of the LGM, the area of the TWFF was covered with snow 
and ice, with accompanying cool temperatures (CLIMAP Project Members, 
1976; Mix et al., 2011). Following the LGM, the retreat of snow and ice, 
together with an increase in temperature, led to a warmer and drier climate that 
likely facilitated the colonization and maintenance of a greater diversity of 
species (Currie et al., 2004). Some species in the region, such as Prunus 
sogdiana, are likely to increase in distribution in wetter conditions (Li et al., 
2011) and may have led to higher levels of diversity at sites where levels of 
precipitation have increased since the LGM. 


4.4 Relationship between tree and herb species diversity 


A recent review proposed that, in temperate forests, the contribution of herbs 
to ecosystem function is context-dependent (Landuyt et al., 2019), while 
previous studies have shown that the relationship between tree and herb 
species diversity may be positive or negative, or non-existent (Berger and 
Puettmann, 2000; Ewald and Freising-Weihenstephan, 2002; Gilliam, 2007; 
Houle, 2007; Molder et al., 2008; Vockenhuber et al., 2011). In the TWFF, 
diversity of tree species was positively associated with that of herb species 
(Figs. 5a—c), supporting previous studies of temperate and deciduous forests 
(Gilliam, 2007; Molder et al., 2008; Vockenhuber et al., 2011), and it was 
unsurprising that we also found a positive relationship between tree and forb 
species diversity (Figs. 5d—f), given the large proportion of forb herbs in this 
area. 

This positive relationship between tree and herb species diversity may be 
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driven by a number of mechanisms. Firstly, the greater diversity of tree species 
leads to more heterogeneous availability of light and soil nutrients, which 
further increases herb species diversity (Huston, 1994; Barbier et al., 2008; 
Vockenhuber et al., 2011). In the TWFF, high levels of herb species diversity at 
the Huocheng study site was mainly attributed to the presence of shade- 
tolerant herbs, such as Impatiens brachycentra and Viola collina (Wang et al., 
2016). Low levels of herb species diversity at the Xinyuan and Gongliu study 
sites were caused by high densities of species with rapid growth rates, such as 
Urtica dioica and Origanum vulgare (Li et al., 2011), that then increase 
competition among slower growing species for light and reduce overall herb 
species diversity. Secondly, nutrient compensation between trees and herbs is 
a key driver of the positive relationship between tree and herb species 
diversity. For example, diversity of herb species in forests of West Virginia is 
regulated by soil calcium (Ca), magnesium (Mg), and potassium (K), while tree 
species diversity is driven by soil phosphorus (P) (Gilliam and Robert, 2003; 
Gilliam, 2007), as indicated by greater concentrations of P in herb leaves than 
in trees and greater concentrations of Ca** in tree leaves than in herbs (Muller, 
2003; Gilliam, 2007). Thus, decomposed P from herbs becomes available for 
tree growth and maintains high levels of species diversity, while releases of 
Ca? from tree leaf litter provide resource for herbs. 


5 Conclusions 


Understanding the diversity relationships between tree and herb species has 
an important implication in forest diversity conservation. In the TWFF, tree 
species diversity was positively associated with herb species diversity due to 
they responded similarly along the latitudinal, elevational and climatic 
gradients, indicating that both groups should be protected to ensure adequate 
conservation of this unique forest ecosystem. However, the underlying 
mechanisms of this positive relationship between tree and herb species 
diversity are unclear, so we recommend that controlled experiments are further 
needed to explore resource competition and nutrient return among and 
between the two species groups. 
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